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Final  Report  on  Seedling  Project  “Suppression  of  Laser  Shot  Noise  Using  Laser- 
Cooled  Optomechanical  Systems”  (Contract  #:  W911NF-09-1-0015)  April  22,  2010 

Jack  Harris,  Department  of  Physics,  Yale  University,  New  Haven,  CT. 

Executive  Summary:  The  goal  of  this  effort  was  to  create  squeezed  light  with 
intensity  fluctuations  15  dB  below  the  shot  noise  level.  This  goal  was  to  be  achieved  by 
integrating  ultrasensitive  MEMS  inside  a  high  finesse  optical  cavity  at  cryogenic 
temperatures.  While  we  have  not  yet  demonstrated  squeezing,  we  have  overcome  what 
proved  to  be  the  most  challenging  technical  hurdle:  coupling  laser  light  into  a  high- 
finesse  cavity  containing  an  ultrasensitive  MEMS  at  T  =  5  K.  Since  this  was  the  major 
technical  novelty  in  this  project,  we  expect  that  this  device  will  be  able  to  demonstrate 
squeezing  in  a  fairly  short  time. 

Background:  The  goal  of  this  effort  was  to  create  laser  light  with  intensity 
fluctuations  below  the  shot  noise  level.  This  type  of  light  (known  as  “squeezed”  light) 
can  improve  the  performance  of  optical  instruments  which  are  presently  limited  by 
photon  shot  noise.  Specific  theoretical  proposals  have  considered  applications  for 
squeezed  light  in  laser  gyroscopes,  trace  gas  detectors,  interferometers,  and 
magnetometers. 

To  date,  squeezed  light  has  been  produced  only  in  proof-of-principle  experiments, 
and  has  not  yet  been  applied  to  any  instruments.  The  present  world  record  for  shot  noise 
suppression  is  10  dB,  but  achieving  this  typically  requires  a  full  optical  table  of 
equipment  and  constant  fine  tuning.  The  principle  challenge  in  producing  squeezed  light 
has  been  meeting  two  simultaneous  requirements  in  a  single  device:  strong  optical 
nonlinearity  (to  alter  the  photon  statistics)  and  low  optical  absorption  (to  avoid  adding 
quantum  noise  to  the  laser  beam).  In  most  nonlinear  optical  materials  these  two  quantities 
are  not  independent,  and  typically  one  must  strike  a  compromise  between  them. 

The  goal  of  our  work  has  been  to  use  a  novel  form  of  optical  nonlinearity  to 
circumvent  this  compromise  and  produce  strongly  squeezed  light  (15  dB  suppression  of 
shot  noise).  Specifically,  we  have  developed  an  ultrasensitive  MEMS  device  which  is 
coupled  via  radiation  pressure  to  an  optical  cavity.  This  optomechanical  device  is 
expected  to  achieve  optical  nonlinearity  because  photons  in  the  cavity  displace  the 
membrane  (via  radiation  pressure),  and  the  displacement  of  the  membrane  in  turn  detunes 
the  cavity. 

The  resulting  intensity-dependent  cavity  detuning  is  equivalent  to  a  Kerr-type 
nonlinearity.  However  this  optomechanical  device’s  nonlinearity  is  achieved  without 
passing  the  light  through  a  solid  medium.  As  a  result,  its  optical  absorption  can  be  much 
lower  than  in  traditional  nonlinear  devices.  This  allows  the  cavity  finesse  to  be  increased, 
leading  to  stronger  optical  nonlinearity. 

The  squeezing  performance  of  such  a  device  has  been  calculated  analytically.  These 
calculations  indicate  that  the  device’s  performance  is  optimized  when  the  cavity  finesse  F 
is  high,  the  membrane’s  mechanical  quality  factor  Q  is  high,  and  the  membrane’s 
Brownian  motion  is  small. 

Accomplishments  During  This  Seedling  Project:  To  meet  these  technical 
requirements,  we  began  this  seedling  project  by  designing  and  constructing  a  high-finesse 


optical  cavity  that  incorporates  a  micromechanical  membrane  and  is  capable  of  operating 
inside  a  cryostat. 

Figure  1  shows  the  basic  schematic  of  the  apparatus.  The  cavity  is  assembled  at 
room  temperature  and  laser  light  is  coupled  into  the  cavity  using  free  space  optics.  Upon 
cooling  from  room  temperature  to  the  cryostat’s  base  temperature,  thermal  contractions 
inevitably  cause  the  cavity  to  become  misaligned.  Provided  this  misalignment  is  small,  it 
can  be  compensated  for  by  adjusting  the  room  temperature  optics  on  the  top  of  the 
cryostat.  The  membrane’s  alignment  to  the  cavity  axis  must  also  be  maintained  during  the 
cool  down.  This  is  accomplished  through  the  use  of  in  situ  cryogenic  translation  stages. 

Since  the  operation  of  a  high-finesse  cavity  inside  a  cryostat  is  a  highly  non¬ 
standard  requirement,  at  the  start  of  this  project  (January  2009  -  March  2009)  we 
conducted  several  incremental  tests  of  this  apparatus’  performance.  First  we  used  our 
cryostat  to  cool  a  high-finesse  cavity  without  a  micromechanical  membrane.  The  purpose 
of  this  was  to  test  whether  we  could  use  free  space  optics  to  maintain  coupling  to  the 
optical  cavity  while  it  cooled.  This  test  was  successful;  thermal  contractions  inside  the 
cryostat  could  be  compensated  for  using  the  room  temperature  optics.  Figure  2A  shows 
measurements  of  this  “empty”  cavity’s  ringdown  at  T  =  4  K.  The  cavity  finesse 
determined  from  this  data  is  in  good  agreement  with  the  value  measured  at  room 
temperature,  indicating  that  the  cavity  mirrors  are  essentially  unaffected  by  operation  at 
low  temperatures.  Similar  results  were  achieved  with  the  fridge  at  T  =  0.4  K,  and  we 
were  able  to  determine  that  laser  powers  up  to  0.1  mW  did  not  heat  the  fridge  by  more 
than  a  few  mK. 

The  second  incremental  test  we  performed  was  to  measure  the  membrane’s  optical 
absorption.  For  simplicity,  these  measurements  were  performed  at  room  temperature.  Our 
previous  work  had  used  non-stochiometric  SiNx  membranes  whose  absorption 
corresponded  to  an  imaginary  index  of  refraction  Im (n)  =  2x1  O'4  (this  data  is  shown  on 
the  left-hand  side  of  Figure  2B).  This  absorption,  while  fairly  low,  meant  that  if  the 
membrane  was  placed  anywhere  but  at  a  node  if  the  intracavity  standing  wave,  the  need 
to  avoid  heating  the  membrane  via  optical  absorption  limited  the  useable  cavity  finesse  to 
-10,000. 

We  had  been  told  by  colleagues  at  Caltech  (O.  Painter,  C.  Regal)  that  stochiometric 
Si3N4  membranes  might  have  much  lower  absorption,  so  during  this  seedling  project  we 
repeated  these  measurements  with  a  stochiometric  membrane,  producing  the  results  on 
the  right-hand  side  of  Figure  2B.  The  cavity’s  finesse  is  modulated  by  only  a  few  percent 
as  the  membrane  is  moved  from  node  to  antinode,  but  since  this  modulation  is  not 
sinusoidal,  it  does  not  appear  to  be  due  to  absorption  in  the  membrane  (which  should  be 
proportional  to  the  overlap  between  the  optical  standing  wave  and  the  membrane).  The 
observed  modulation  of  the  cavity  finesse  is  believed  to  be  due  to  the  scattering  of  light 
between  different  transverse  modes.  Although  such  scattering  can  limit  the  usable 
finesse,  it  does  not  result  in  heating  of  the  membrane.  Based  on  the  data  in  Figure  2B  we 
estimate  that  Im (n)  <  2xl0"4  in  these  stochiometric  membranes.  These  stochiometric 
membranes  achieve  the  same  high  mechanical  quality  factor  ( Q  =  lxlO6)  as  the 
nonstochiometric  membranes.  These  results  are  incorporated  in  a  paper  which  has  just 
been  accepted  for  publication  in  Nature  Physics  (currently  available  as 
ArXiv:  1002.4158). 


With  these  preliminary  tests  conducted,  we  proceeded  to  design  and  build  the  actual 
device  (April  2009  -  May  2009).  The  cryogenic  cavity  demonstrated  in  Figure  2 A  was 
modified  to  include  a  membrane,  as  well  as  cryogenic  translation  stages  for  aligning  the 
membrane  within  the  cavity.  The  final  apparatus  is  shown  in  Figure  3,  which  gives  both 
an  external  view  of  the  room  temperature  optics  and  an  internal  view  of  the  device  itself. 

In  cooling  down  the  modified  cavity  (which  now  incorporated  the  membrane)  we 
found  that  thermal  contractions  (or  other  effects)  left  the  cavity  too  misaligned  to  be 
corrected  using  the  room  temperature  optics.  This  does  not  require  a  very  large 
misalignment  -  once  the  cavity  axis  has  shifted  enough  that  it  no  longer  points  out  the  1 
cm  diameter,  2  m  long  clear  shot  tube,  it  is  no  longer  possible  for  the  room  temperature 
optics  to  couple  a  laser  beam  to  the  cavity. 

The  period  June  2009  -  March  2010  was  spent  debugging  this  problem.  The 
solution  involved  a  careful  redesign  of  the  cryogenic  optical  setup  to  minimize 
asymmetry,  strains,  and  dissimilar  materials.  In  early  March,  we  were  able  to  cool  down 
our  device  while  maintaining  full  optical  coupling  to  the  cavity  and  full  control  over  the 
membrane’s  position  inside  the  cavity. 

Figure  4  shows  a  measurement  of  the  cavity’s  optical  ringdown  (taken  at  T  =  5  K) 
with  the  membrane  at  a  node  of  the  cavity  mode  (upper  panel)  and  at  an  antinode  (lower 
panel).  As  expected,  the  measured  cavity  finesse  varies  by  ~1%,  indicating  that  the 
membrane  is  not  absorbing  an  appreciable  amount  of  the  intracavity  light. 

Figure  5  shows  a  measurement  of  the  membrane’s  mechanical  resonance  near  523 
kHz.  The  mechanical  Q  factor  appears  low  because  this  data  was  taken  with  helium 
exchange  gas  inside  the  fridge’s  vacuum  space.  We  have  since  removed  the  exchange  gas 
and  have  begun  characterizing  the  membrane  in  vacuum. 

The  most  recent  data  we  have  taken  is  shown  in  Figure  6,  which  plots  the 
fluctuations  in  the  intensity  of  a  laser  beam  reflecting  off  of  the  cryogenic 
optomechanical  device  described  above.  The  goal  of  this  project  was  to  have  the 
fluctuations  of  this  beam  suppressed  below  the  shot  noise  level.  As  described  above,  we 
have  not  yet  achieved  this  (mostly  owing  to  the  delay  in  debugging  the  cryogenic 
operation  of  the  experiment),  but  with  the  apparatus  working  in  the  last  month,  we  are 
rapidly  approaching  this  goal. 

In  Figure  6,  the  black  curve  (lowermost)  is  the  detector  noise  measured  in  the 
absence  of  a  laser  beam.  The  dashed  curve  (just  above  the  detector  noise)  is  the  expected 
shot  noise  for  the  laser  power  used  in  the  remaining  data  traces.  As  can  be  seen,  it  is 
above  the  detector  noise  and  so  should  be  observable.  The  blue  trace  in  Figure  6 
corresponds  to  a  laser  beam  which  is  deliberately  detuned  from  the  cavity  resonance,  and 
so  is  simply  reflected  from  the  front  mirror  of  the  cavity  without  interacting  with  either 
the  cavity  or  the  membrane.  This  serves  a  baseline  measurement  of  our  laser  noise.  As 
can  be  seen,  it  is  roughly  an  order  of  magnitude  above  the  expected  shot  noise. 
Presumably  this  represents  technical  noise  in  our  laser;  in  future  measurements  this  will 
be  passively  filtered  using  a  reference  cavity,  which  is  already  assembled  and  tested  in 
our  lab. 

The  green,  gray,  and  red  traces  correspond  to  data  taken  with  the  laser  locked  to  the 
cavity  resonance  (and  so  interacting  with  the  cavity  and  the  membrane).  In  the  gray 
trace,  the  membrane  is  positioned  at  a  node  of  the  intracavity  field;  as  a  result,  we  expect 
it  to  interact  minimally  with  the  laser  beam  and  so  to  produce  very  little  optomechanical 


effect.  In  the  green  and  red  traces  the  membrane  is  displaced  slightly  from  the  node,  into 
a  region  where  the  optical  standing  wave  overlaps  more  strongly  with  the  membrane.  At 
present  we  have  not  observed  any  impact  of  the  membrane’s  motion  on  the  laser  noise, 
but  we  are  presently  working  on  this.  One  important  issue  is  that  the  data  in  Figure  6  was 
taken  with  the  laser  tuned  precisely  to  the  cavity  resonance,  whereas  the  strongest 
optomechanical  effects  are  expected  when  the  laser  is  detuned  by  roughly  one  linewidth. 

It  should  be  possible  to  modify  the  laser  detuning  and  membrane  position  in  situ 
with  the  present  set  up.  Once  this  is  accomplished  and  other  minor  technical  issues  are 
debugged,  we  expect  to  see  the  membrane’s  response  suppress  the  laser  noise.  We  expect 
that  this  will  require  a  few  more  months,  and  that  the  single  greatest  technical  challenge 
(maintaining  optical  coupling  to  the  cryogenic  cavity)  has  been  overcome  during  the 
period  of  this  seedling. 

We  note  that  in  parallel  with  these  experimental  efforts  we  have  collaborated  with 
theory  colleagues  at  Yale  (Steven  Girvin  and  his  group)  to  consider  the  possibilities  for 
generating  squeezed  light  via  different  types  of  optomechanical  coupling.  We  have 
developed  a  theory  for  the  squeezing  produced  by  a  quadratic  optomechanical  coupling, 
which  our  devices  are  capable  of  realizing  (in  addition  to  the  usual  linear  coupling).  This 
work  has  been  submitted  to  Physical  Review  Letters. 

Publications  Resulting  From  This  Work: 

oNunnenkamp,  K.  Borkje,  J.  G.  E.  Harris,  S.  M.  Girvin,  Cooling  and  squeezing  via 
quadratic  optomechanical  coupling,  ArXiv:  1004.25 10  (submitted  to  Physical 
Review  Letters,  2010). 

o  J.  C.  Sankey,  C.  Yang,  B.  M.  Zwickl,  A.  E.  Jayich,  J.  G.  E.  Harris,  Strong  and 
tunable  nonlinear  optomechanical  coupling  in  a  low-loss  system, 
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FIGURE  1:  Schematic  of  a  cryogenic  high-finesse 
cavity  with  an  integrated  micromechanical 
membrane.  Laser  light  is  coupled  into  the  cavity 
from  room  temperature  via  free  space  optics. 


stage 


FIGURE  2:  Preliminary  characterization 


Figure  2A:  Optical  ringdown  of  an  empty 
cavity  (i.e.,  without  a  membrane)  inside 
our  cryostat  at  T  =  4K.  The  fit  is  to  a 
single  exponential  and  returns  a  value  for 
the  finesse  of  71 ,000,  equivalent  to  the 
value  measured  at  room  temperature. 
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Figure  2b:  Comparison  of  absorption  for  a  nonstochiometric  membrane  (left  panel)  and  a  stochiometric 
membrane  (left  panel).  The  nonstochiometric  membrane  absorbs  strongly  whenever  it  is  away  from  a  node  of 
the  intracavity  standing  wave  (arrows),  resulting  in  a  strong  decrease  of  the  cavity  finesse.  The  dashed  line 
shows  the  cavity  finesse  when  the  membrane  is  removed.  The  solid  line  is  a  fit  to  the  data  and  gives  Im (n)  =  2 
xIO'4.  Stochiometric  membranes,  by  contrast,  absorb  very  little  light  (right  panel).  For  this  measurement  the 
empty  cavity  finesse  was  50,000.  The  small  modulations  in  the  finesse  reflect  scattering  between  transverse 
modes,  not  absorption. 


Figure  3:  The  experimental  cryostat.  Left:  photograph  of  the  exterior  of  the  cryostat  (blue).  An  optical 
breadboard  (black)  is  mounted  to  the  top  of  the  cryostat,  allowing  laser  light  to  be  aligned  down  the  cryostat’s 
clear-shot  tubes.  The  laser  light  is  delivered  to  the  cryostat  via  optical  fibers  (yellow).  Right:  photograph  of  the 
interior  of  the  cryostat,  showing  the  3He  cold  plate  (gold),  the  cavity,  and  the  piezoelectric  translation  stages 
which  allow  the  membrane  to  be  positioned  within  the  cavity  in  situ. 


Figure  4:  Optical  ringdown  of  the  cavity  containing 
a  membrane  at  T  =  5  K.  Upper  panel:  ringdown 
with  the  membrane  at  a  node  of  the  optical 
standing  wave.  The  red  curve  is  a  fit  to  a  single 
exponential  giving  a  cavity  finesse  80,000.  Lower 
panel:  ringdown  with  the  membrane  at  an  antinode 
of  the  optical  standing  wave.  The  red  curve  is  a  fit 
to  a  single  exponential  giving  a  cavity  finesse 
80,000. 


Figure  5:  Measurement  of  the  mechanical 
resonance  of  the  membrane  inside  the  cavity  at  T 
=  5  K.  The  low  quality  factor  is  due  to  the  fact  that 
this  data  was  taken  with  helium  exchange  gas 
inside  the  fridge. 
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Figure  6:  Preliminary  measurement  of  laser  noise  reflected  from  a  cryogenic  (T  =  5  K)  optomechanical 
cavity.  The  expected  shot  noise  (dashed  line)  is  well  above  the  detector  dark  noise  (black).  However  it  is 
at  present  obscured  by  technical  noise  from  the  laser  (blue).  Tuning  the  laser  onto  resonance  with  the 
optomechanical  cavity  (gray,  green,  and  red)  produces  small  changes  in  the  observed  noise,  but 
substantial  optimization  of  the  input  laser  noise  as  well  as  laser  power  and  detuning  relative  to  the  cavity 
remains  to  be  done  in  the  near  term. 
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Abstract 

A  major  goal  in  optomechanics  is  to  observe  and  control  quantum  behavior  in  a  system 
consisting  of  a  mechanical  resonator  coupled  to  an  optical  cavity.  Work  towards  this  goal  has 
focused  on  increasing  the  strength  of  the  coupling  between  the  mechanical  and  optical  degrees 
of  freedom;  however,  the  form  of  this  coupling  is  crucial  in  determining  which  phenomena  can 
be  observed  in  such  a  system.  Here  we  demonstrate  that  avoided  crossings  in  the  spectrum  of  an 
optical  cavity  containing  a  flexible  dielectric  membrane  allow  us  to  realize  several  different  forms 
of  the  optomechanical  coupling.  These  include  cavity  detunings  that  are  (to  lowest  order)  linear, 
quadratic,  or  quartic  in  the  membrane’s  displacement,  and  a  cavity  finesse  that  is  linear  in  (or 
independent  of)  the  membrane’s  displacement.  All  these  couplings  are  realized  in  a  single  device 
with  extremely  low  optical  loss  and  can  be  tuned  over  a  wide  range  in  situ ;  in  particular,  we  find 
that  the  quadratic  coupling  can  be  increased  three  orders  of  magnitude  beyond  previous  devices. 

As  a  result  of  these  advances,  the  device  presented  here  should  be  capable  of  demonstrating  the 
quantization  of  the  membrane’s  mechanical  energy. 

Nearly  all  optomechanical  systems  realized  to  date  can  be  characterized  by  a  linear  relationship 
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between  the  optical  cavity’s  detuning  uj{x)  and  the  displacement  of  the  mechanical  element  x.1 
In  the  classical  regime  this  “linear”  optomechanical  coupling  has  enabled  powerful  laser  cooling 
and  sensitive  displacement  readout  of  the  mechanical  element.2-7  As  uJ  =  duo/dx  increases  this 
linear  coupling  becomes  stronger,  and  it  should  become  possible  to  observe  quantum  effects  such 
as  laser-cooling  to  the  mechanical  ground  state,8,9  quantum-limited  measurements  of  force  and 
displacement,10, 11  and  the  production  of  squeezed  light.12  In  the  quantum  regime,  however,  the 
form  of  the  optomechanical  coupling  plays  a  crucial  role  in  determining  which  phenomena  are 
observable.  For  example,  linear  coupling  provides  a  continuous  readout  of  x,  and  so  precludes  a 
direct  measurement  of  one  of  the  most  striking  features  associated  with  the  quantum  regime:  the 
quantization  of  the  mechanical  oscillator’s  energy. 

One  device  that  has  demonstrated  a  nonlinear  optomechanical  coupling  consists  of  a  thin  dielec¬ 
tric  membrane  placed  inside  a  Fabry-Perot  cavity.13  With  the  membrane  positioned  at  a  node  (or 
antinode)  of  the  intracavity  standing  wave,  w(x)  oc  x1  to  lowest  order.  This  “quadratic”  optome¬ 
chanical  coupling  is  compatible  with  a  quantum  nondemolition  (QND)  readout  of  the  membrane’s 
energy  Hm  —  where  ujm  is  the  membrane’s  resonant  frequency  and  nm  is  the  membrane’s 

phonon  number.  Two  distinct  schemes  have  been  proposed  for  using  this  quadratic  coupling  to 
demonstrate  the  quantization  of  the  membrane’s  energy.  Both  schemes  assume  the  membrane  is 
laser-cooled  to  mean  phonon  number  (nm)  <  1.  Numerical  estimates  indicate  this  level  of  cooling 
should  be  feasible  for  the  device  described  here,  provided  it  is  pre-cooled  cryogenically.8,9 

The  goal  of  the  first  scheme  is  to  monitor  nm  with  resolution  sufficient  to  observe  individual 
quantum  jumps.  This  is  not  feasible  with  the  devices  demonstrated  to  date,  and  would  require 
substantial  improvements  to  the  quadratic  coupling  strength  cj"  =  d2u)/dx 2,  the  membrane’s  optical 
absorption,  and  the  membrane’s  mechanical  properties.13 

In  the  second  scheme,  the  laser-cooled  membrane  would  be  mechanically  driven  from  the  ground 
state  to  a  large-amplitude  coherent  state  with  (nm)  ^$>  1.  The  quadratic  coupling  would  then  be 
used  to  monitor  the  membrane’s  energy  with  resolution  sufficient  to  resolve  fluctuations  oc  yj  (nm) 
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corresponding  to  the  shot  noise  of  the  membrane’s  phonons.  Detailed  calculations14  show  that  this 
second  scheme  is  considerably  less  demanding  than  the  first,  though  it  still  requires  substantial 
improvements  to  uj"  and  the  membrane’s  optical  absorption. 

Here  we  demonstrate  an  optomechanical  device  in  which  uj"  is  increased  by  at  least  three  orders 
of  magnitude,  while  the  membrane’s  optical  absorption  is  substantially  decreased.  The  device 
satisfies  the  requirements  for  observing  phonon  shot  noise  at  a  bath  temperature  T  —  300  mK,  and 
represents  substantial  progress  toward  observing  individual  quantum  jumps.  The  improvement 
in  uj"  is  achieved  by  exploiting  the  full  spectrum  of  the  optical  cavity’s  transverse  modes,  which 
exhibits  numerous  avoided  crossings  as  a  function  of  the  membrane’s  position.  These  crossings  were 
not  considered  in  previous  work,  which  assumed  a  one-dimensional  model  of  the  cavity  and  only  a 
single  transverse  optical  mode.13,15 

In  addition  to  increased  uj"  ,  we  demonstrate  considerable  flexibility  within  a  single  device:  (i) 
uj"  can  be  varied  in  situ  by  adjusting  the  position  and  tilt  of  the  membrane;  (ii)  it  is  possible  to 
tune  uj"  to  zero,  thereby  realizing  a  purely  quartic  optomechanical  coupling  uj{x)  oc  x4;  (iii)  the 
gradient  of  the  cavity  relaxation  n!  =  dnjdx  can  be  tuned  over  a  wide  range  or  set  to  zero;  and 
(iv)  cavity  modes  with  different  forms  of  uj[x)  (e.g.,  linear  and  quadratic)  can  be  simultaneously 
addressed  using  multiple  laser  frequencies,  allowing  for  simultaneous  laser  cooling  and  QND  energy 
readout.  Each  type  of  coupling  offers  distinct  functionality,  and  together  they  represent  a  new  set 
of  tools  for  observing  and  controlling  quantum  effects  in  optomechanical  systems.  We  find  that  the 
features  in  the  cavity  spectrum  responsible  for  these  couplings  are  reproduced  by  a  straightforward 
theoretical  model,  allowing  for  optimization  of  future  devices. 

1  Strong  purely-quadratic  optomechanical  coupling 

Our  optomechanical  system  is  shown  schematically  in  Fig.  la.  Two  fixed  end  mirrors  (radius 
of  curvature  R  —  5  cm)  form  a  Fabry-Perot  cavity  with  free  spectral  range  2.374  GHz  (length 
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L  —  6.313  cm)  and  empty-cavity  finesse  F  —  50,  000.  A  1-mm-square  flexible  Si3N4  membrane 
of  nominal  thickness  t  —  50  nm  and  real  refractive  index  R e[n\  =  2.0  is  placed  near  the  cavity 
waist.  The  membrane  is  mounted  on  a  motorized  stage,  providing  control  over  the  membrane’s 
coarse  position  along  the  cavity  axis  (x)  as  well  as  its  tilt  about  the  two  transverse  axes  (y  and  z). 
Piezoelectric  transducers  allow  for  nanometer-scale  displacements  along  x. 

Figure  lb  shows  the  cavity’s  transmission  spectrum  as  a  function  of  membrane  position  x.  The 
cosine-like  detuning  curves  are  similar  to  those  demonstrated  previously,13  and  achieve  ou" /27t  =  30 
kHz/nm2  at  their  extrema.  The  data  in  Fig.  lb  were  taken  with  the  laser  coupled  to  several  of  the 
cavity’s  lower-order  transverse  modes,  and  a  number  of  apparent  crossings  between  these  modes 
can  be  seen.  We  focus  specifically  on  the  region  highlighted  by  the  dotted  box  in  Fig.  lb,  where 
the  TEMqo  “singlet”  and  the  TEM{2o,  11,02}  “triplet”  modes  cross.  Figure  lc  shows  this  region  in 
greater  detail  with  the  membrane  (i)  lying  in  the  y-z  plane  and  (ii)  tilted  by  0.4  mrad  about  the 
z-axis.  The  behavior  in  Fig.  lc(ii)  is  ubiquitous  among  multiplets  of  higher-order  transverse  modes: 
tilting  the  membrane  lifts  the  multiplet  degeneracy,  with  modes  of  greater  spatial  extent  along  the 
membrane  slope  (y  in  this  case)  perturbed  the  most. 

The  central  result  of  this  article  is  illustrated  in  Fig.  Id,  which  shows  a  high-resolution  scan 
of  the  region  indicated  by  the  dotted  box  in  Fig.  lc.  This  data  demonstrates  that  the  apparent 
crossings  between  the  various  optical  modes  are  avoided,  and  that  at  their  anticrossings  oj(x)  is 
purely  quadratic  with  oj"/2tt  substantially  greater  than  30  kHz/nm2.  For  example,  the  TEM20 
-  TEMqo  crossing  in  Fig.  Id  shows  (jj"/2tt  =  4.5  MHz/nm2  (dashed  white  lines).  The  TEM02  - 
TEMqo  crossing  is  not  resolved  in  the  main  body  of  Fig.  Id,  but  the  line-scan  in  the  inset  shows  a 
splitting  ~  200  kHz  between  these  modes.  Assuming  the  crossing  is  of  the  usual  hyperbolic  form 
u(x)  =  ±y/ ( (jj'x )2  +  where  uJ  is  the  slope  of  u(x)  far  from  the  crossing  and  2us  is  the  gap  at 
the  degeneracy  point,  the  data  in  the  inset  indicate  /2n  >  30  MHz/nm2,  which  is  three  orders 
of  magnitude  greater  than  previously  demonstrated.13 

We  note  that  if  the  membrane  is  positioned  so  that  two  modes  realize  a  quadratic  coupling, 
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other  modes  will  still  realize  a  linear  coupling.  In  Fig.  Id  this  occurs  at  the  position  x  =  0  nm: 
here  the  eigenmodes  formed  by  the  TEM02  and  TEMqo  modes  exhibit  quadratic  coupling  while 
the  TEMn  mode’s  coupling  is  linear.  This  means  that  lasers  tuned  to  different  eigenmodes  can 
simultaneously  exploit  different  forms  of  the  optomechanical  coupling. 

We  can  understand  the  origin  of  these  features  by  noting  that  avoided  crossings  generally 
reflect  a  broken  symmetry  that  prevents  eigenmodes  from  becoming  degenerate.  An  ideal  empty 
Fabry-Perot  cavity  possesses  symmetry  that  allows  degeneracy  between  transverse  modes,  but  in 
our  device  we  expect  this  symmetry  to  be  broken  for  two  reasons:  the  curved  wavefronts  of  the 
cavity  modes  may  not  overlap  perfectly  with  the  flat  membrane  (e.g.,  if  the  membrane  is  tilted 
or  displaced  from  the  cavity  waist),  and  the  empty  cavity  itself  may  be  slightly  asymmetric  (e.g., 
owing  to  imperfect  form  of  the  end  mirrors). 

To  make  a  quantitative  analysis  of  the  cavity  spectrum  and  the  features  in  Fig.  lb-d,  we  devel¬ 
oped  a  perturbative  solution  of  the  Helmholtz  equation  to  calculate  the  eigenmodes  and  eigenfre- 
quencies  of  a  symmetric  optical  cavity  into  which  a  dielectric  slab  is  placed  at  an  arbitrary  location 
and  tilt.  As  described  elsewhere,16  the  empty-cavity  eigenfrequencies  are  perturbed  by  an  amount 
proportional  to  the  eigenvalues  of  the  matrix  V,  where  Vij  oc  fff^i(x,y,z)^j(x,y,z)dxdydz. 
Here  V’fc  is  the  kth  unperturbed  eigenmode  of  the  cavity  and  the  integral  is  taken  over  the  volume 
of  the  membrane. 

Applying  this  theory  to  the  four  cavity  modes  of  interest  (the  singlet  and  triplet  modes)  quan¬ 
titatively  reproduces  the  large-scale  (^GHz)  sinusoidal  shape  of  u(x)  seen  in  Fig.  lb  if  we  assume 
Re[n]=  2.0  and  t  —  39  nm  (black  lines).  The  discrepancy  between  the  fitted  and  nominal  values  of  t 
presumably  reflects  a  combination  of  fabrication  tolerances  (^  5  nm)  and  the  limits  of  a  first-order 
perturbation  theory  that  includes  only  four  eigenmodes.  Nonetheless,  the  intermediate-scale  (^10 
MHz)  features  in  Fig.  lc(i-ii),  such  as  the  lifting  of  the  triplet’s  degeneracy,  are  also  quantitatively 
reproduced  (Fig.  lc(iii-iv)).  Small-scale  features  such  as  avoided  crossings  agree  with  the  model 
reasonably  well  and  are  discussed  below. 
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2  Tunability  of  quadratic  coupling 


The  second  result  of  this  article  is  that  u"  can  be  tuned  over  a  wide  range  by  moving  the  membrane 
along  the  x  axis.  This  tunability  is  important  because  in  some  situations  it  may  be  desirable  to 
decrease  uo"  in  order  to  relax  other  experimental  constraints.  For  example,  if  uos  (oc  l/cc/7)  is  small 
enough  the  membrane’s  motion  may  result  in  non-adiabatic  transfer  of  light  between  the  two  cavity 
modes  via  Landau-Zener-Stiickelberg-like  transitions.17 

The  tunability  of  u"  is  illustrated  in  Fig.  2a-b,  which  each  show  six  avoided  crossings  between 
the  singlet  and  triplet  modes.  When  the  membrane  is  at  the  cavity  waist  (Fig.  2a)  the  upper  gaps 
(triangles)  are  open  and  the  lower  gaps  (squares)  are  closed.  When  the  membrane  is  displaced  500 
pm  (Fig.  2b)  from  the  waist,  the  two  lower  gaps  open,  the  upper  right  gap  opens  further,  and  the 
upper  left  gap  closes.  The  full  dependence  of  uos  on  membrane  position  is  shown  in  Fig.  2c. 

The  perturbative  model  (Fig.  2c  inset)  reproduces  the  linear  dependence  of  uos{x)  as  well  as  the 
slope  dujs/dx  measured  for  each  of  the  six  gaps.  The  model  differs  from  the  data  by  a  constant 
offset  ^3  —  4  MHz  for  the  middle  and  upper  gaps  (triangles  and  crosses  in  Fig.  2c),  which  we 
attribute  to  asymmetry  in  the  cavity  that  is  not  included  in  the  model.  We  find  below  (Fig.  3c) 
that  it  is  possible  to  compensate  for  this  intrinsic  asymmetry  by  adjusting  the  axis  of  the  membrane 
tilt. 

3  Tunable  coupling  between  motion  and  optical  relaxation 

In  addition  to  creating  large  uo" ,  mixing  between  the  cavity  modes  also  causes  the  cavity  loss  rate  ft  to 
vary  with  x.  Figure  3a(i)  shows  k{x)  for  the  singlet  and  triplet  modes  (the  corresponding  detunings 
are  shown  in  Fig.  3a(ii)).  Far  from  the  crossings  each  cavity  mode  has  a  different  ft,  reflecting  the 
different  overlap  between  each  mode’s  transverse  profile  and  the  mirrors’  inhomogeneities.  At  each 
avoided  crossing  the  two  modes  swap  their  value  of  ft,  leading  to  a  large  ft'  =  dhi/dx.  The  dashed 
line  in  Fig.  3a  indicates  a  crossing  at  which  hi'  >  600  kHz/nm. 
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Note  that  this  data  was  taken  with  the  membrane’s  tilt  axis  rotated  by  45°  to  (y  +  z)/y/2,  and 
the  transverse  eigenmodes  have  followed  (inset  camera  images).  The  ability  to  rotate  the  transverse 
mode  profiles  provides  some  control  over  which  portions  of  the  mirrors  the  optical  modes  sample, 
and  enables  us  to  tune  n!  in  situ .  In  Fig.  3b  we  have  rotated  the  membrane  tilt  axis  back  to  5, 
with  the  result  that  the  same  crossing  shows  kI  =  — 10  kHz/nm.  Significantly,  the  sign  of  n!  has 
changed,  implying  that  at  an  intermediate  tilt  axis  n!  —  0. 

Gradients  in  k  can  have  several  consequences.  A  linear  variation  of  k  with  x  could  preclude 
a  QND  measurement  of  nm,  so  the  ability  to  tune  kI  to  zero  is  appealing.  Separately,  it  has 
been  predicted  that  large  values  of  k!  can  be  used  to  laser  cool  the  membrane  to  its  ground  state 
even  when  the  device  is  not  in  the  resolved  sideband  regime.18  However  numerical  estimates  show 
that  even  with  kI  —  600  kHz/nm,  the  usual  optomechanical  coupling8,9  is  more  promising  for 
laser-cooling  the  devices  described  here. 

Variations  in  k  can  also  arise  from  optical  absorption  in  the  membrane.  Previous  measurements 
using  non-stoichiometric  SiNx  membranes  found  that  variations  in  n  were  proportional  to  the 
overlap  of  the  intracavity  standing  wave  and  the  membrane,  and  were  consistent  with  Im(n)  ~ 
2  x  10-4  for  A  =  1064  nm.15  Subsequent  work  found  lower  absorption  in  stoichiometric  Si3N4 
membranes,  with  Im(n)  <  10-5  for  A  =  985  nm.19  Figure  3c  shows  n  for  the  TEMqo  mode 
(using  a  Si3N4  membrane  and  A  =  1064  nm)  as  x  is  varied  over  a  few  A.  The  small-scale  (^C  A) 
variations  in  n(x)  are  reproducible  and  periodic,  and  presumably  arise  from  mixing  with  higher- 
order  modes.  To  estimate  the  contribution  to  n(x)  from  absorption  in  the  membrane  (as  opposed 
to  mixing  with  higher-order  modes),  we  extract  the  Fourier  component  of  k(x)  corresponding  to 
the  overlap  between  the  membrane  and  the  intracavity  standing  wave.  This  sets  an  upper  limit  of 
Im(n)  <  1.5  x  10-6  at  A  =  1064  nm.  The  lower  optical  absorption  observed  here  should  enable  the 
use  of  cavities  with  higher  F  while  decreasing  the  quantum  noise  in  the  cavity  and  the  heating  of 
the  membrane. 
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m 


F  J  '/2n  Um/ 27t  Qm 


Quantum  Jumps: 

300,000 

10  MHz/nm2 

100  kHz 

1.2  x  107 

50  pg 

£(°)  = 

1.1 

Phonon  Shot  Noise: 

50,000 

0.9  MHz/nm2 

1  MHz 

1.2  x  107 

40  ng 

5  = 

7.9 

Table  1:  Comparison  of  parameters  (cavity  length  L,  finesse  F,  quadratic  coupling  c */',  membrane 
frequency  ujm,  quality  factor  Qm ,  and  mass  m)  to  observe  energy  quantization  in  a  mechanical 
resonator  for  the  two  schemes.  Both  cases  assume  5  pW  of  1064-nm  light  incident  on  a  cavity  of 
length  L  —  6.313  cm  in  a  3He  cryostat,  the  membrane  positioned  within  0.5  pm  of  the  avoided 
crossing,  and  the  mechanical  motion  laser-cooled  to  tit  <  0.2  from  a  starting  temperature  T  = 
300  mK.  For  the  phonon  shot  noise  experiment,  the  estimate  assumes  the  membrane  is  driven  to  1 
nm  amplitude  of  motion. 

4  Feasibility  of  observing  mechanical  energy  quantization 

To  estimate  the  feasibility  of  observing  energy  quantization  in  the  membrane,  we  first  use  the  ex¬ 
pressions  derived  elsewhere13  for  E®,  the  signal-to- noise  ratio  for  measurement  of  a  single  quantum 
jump.  Assuming  a  single-port  cavity,20  we  estimate  =  1.1  for  a  device  with  the  parameters 
listed  in  the  first  row  of  Table  1.  This  scheme  requires  higher  finesse  and  smaller  values  of  m  and 
ujm  (which  could  be  realized  by  patterning  the  membrane  into  a  free  standing  pad  supported  by 
narrow  beams),  as  well  as  cryogenic  pre-cooling  to  T  —  300  mK.  The  key  advances  presented  in 
this  paper  regarding  the  realization  of  this  scheme  are  the  demonstration  of  (i)  sufficiently  large  uj" 
and  (ii)  sufficiently  low  optical  loss  to  meet  these  requirements. 

The  requirements  for  observing  phonon  shot  noise  in  the  membrane  are  described  in  Ref.  [14]. 
The  ratio  between  the  phonon  shot  noise  signal  and  the  measurement  imprecision  for  such  a  mea¬ 
surement  is  S  =  8 nmu/rX(°),  where  nr  =  k^T/huJm.  The  parameters  listed  in  the  second  row  of 
Table  1  result  in  S  =  7.9.  Significantly,  these  parameters  correspond  to  the  device  demonstrated 
here;  the  value  of  uo" /27T  =  0.9  MHz/nm2  corresponds  to  the  lower  avoided  crossings  resolved  in  Fig. 
2b,  and  the  amplitude  of  the  membrane’s  motion  (xq  =  1  nm)  is  well  below  the  onset  of  dynamical 


bistability.21  A  quality  factor  Qm  —  1.2  x  107  was  demonstrated  in  similar  membranes  at  T  —  300 
mK.21 

5  Purely-quartic  optomechanical  coupling 

The  final  point  of  this  article  is  to  demonstrate  a  new  type  of  optomechanical  nonlinearity:  quartic 
(x4)  coupling.  Figure  4a  shows  that  when  the  membrane  tilt  is  increased  to  ^  1  mrad,  uo(x)  for  the 
TEM20  mode  undergoes  a  smooth  transition  from  uo"  >  0  (i)  to  uj"  <  0  (iii);  between  these  limits 
c jj "  —  0  (ii)  and  u{x)  oc  x4  to  lowest  order.  A  similar  transition  is  visible  in  the  faint  background 
modes  of  Fig.  2a  and  Fig.  4b  as  a  function  of  mode  index. 

This  form  of  cj(x)  can  be  used  to  realize  an  optomechanical  coupling  described  (in  the  rotating 
wave  approximation)  by  the  Hamiltonian  term  ffcoup  =  /kj(4)x4pfn7n^,  with  =  d4u4/dx 4, 
rrzpf  =  yj h/2m(jjrn ,  and  n7  the  intracavity  photon  number.  This  type  of  coupling  can  be  used,  for 
example,  to  prepare  Schrodinger  cat  states  in  the  membrane.22,23 

While  the  quartic  coupling  in  Fig.  4a  is  quite  weak  /2i r  =  0.4  Hz/nm4),  it  may  be  possible 

to  increase  using  avoided  crossings.  For  example  the  interaction  between  the  triplet  and 
quintuplet  modes  (Fig.  4b)  shows  avoided  crossings  in  which  uj"  changes  sign.  In  analogy  with  the 
tunability  of  oj"  demonstrated  in  Fig.  2c,  we  expect  that  careful  arrangement  of  the  membrane  tilt 
and  position  will  allow  some  of  the  crossings  in  Fig.  4b  to  be  purely  quartic  with  a  substantially 
larger  u^4\ 

6  Summary 

In  summary,  we  have  demonstrated  an  optomechanical  device  in  which  the  strength  and  the  form 
of  the  optomechanical  coupling  can  be  tuned  over  a  wide  range  in  situ.  We  have  demonstrated 
control  over  whether  the  optical  cavity  detuning  is  (to  lowest  order)  linear,  quadratic,  or  quartic 
in  the  displacement  of  a  micromechanical  membrane,  and  shown  that  the  quadratic  coupling  is 
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three  orders  of  magnitude  stronger  than  previously  demonstrated.  This  device  also  demonstrates 
extremely  low  optical  loss,  and  an  optical  loss  gradient  that  can  be  tuned  to  zero.  These  represent 
important  advances  in  the  ongoing  effort  to  observe  and  manipulate  quantum  behavior  in  a  solid 
mechanical  oscillator. 

In  particular,  the  combination  of  low  optical  loss  and  strong  quadratic  coupling  demonstrated 
here  should  enable  the  observation  of  the  membrane’s  energy  quantization  without  further  im¬ 
provements  to  the  present  device.  This  combination  will  also  enable  other  functionalities  related 
to  quadratic  coupling,  including  dispersive  QND  readout  of  the  intracavity  photon  number,  two- 
phonon  laser  cooling,  conditional  squeezing  between  the  reflected  light  and  the  membrane’s  motion, 
and  various  types  of  passive  optical  squeezing.15,24 

7  Methods 

All  measurements  were  performed  at  room  temperature  and  pressure  <  10-5  Torr.  The  end  mirrors 
were  clamped  to  an  invar  spacer,  and  a  mount  equipped  with  three  motorized  actuators  (for  tilting 
and  displacing  the  membrane)  were  mounted  to  this  spacer.  Two  small  piezoelectric  elements 
provided  finer  positioning  of  the  membrane  along  the  cavity  axis.  Spectroscopy  was  performed  by 
sweeping  the  frequency  of  a  continuous  wave  Nd:YAG  laser  from  low  from  low  to  high  and  stepping 
the  membrane  position  with  the  piezo  elements.  The  cavity’s  optical  loss  was  measured  via  cavity 
ringdown. 

8  Author  information 

Correspondence  and  requests  for  materials  should  be  addressed  to  J.G.E.H. 
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Figure  1:  Avoided  crossings  in  cavity  spectra,  a  Schematic  of  our  setup,  b  Cavity  transmission 
as  a  function  of  laser  detuning  and  membrane  displacement.  The  membrane  is  positioned  near 
the  cavity  waist,  and  the  input  laser  is  coupled  to  the  lower-order  transverse  cavity  modes.  The 
transverse  modes  corresponding  to  the  strongest  transmission  peaks  are  labeled.  Images  of  each  of 
these  modes  (captured  with  an  infrared  camera)  are  shown.  Solid  lines  show  model  results.  c(i,ii) 
Refined  scans  of  the  box  labeled  “c”  in  a  with  the  membrane  (i)  aligned  and  (ii)  tilted  by  0.4  mrad. 
(iii)  and  (iv)  show  the  corresponding  predictions  of  the  model,  d  Refined  scan  of  the  box  labeled 
“d”  in  c(ii) . 
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Figure  2:  Tunability  of  avoided  crossings  for  the  membrane  tilted  0.48  mrad  about  5.  a  Transmission 
spectrum  for  the  membrane  positioned  near  the  cavity  waist  {x  —  0  fim).  b  Transmission  spectrum 
for  the  membrane  displaced  by  -0.5  mm.  c  Magnitudes  of  the  gaps  labeled  in  a  &  b  plotted  as 
a  function  of  the  membrane’s  displacement.  Error  bars  indicate  the  frequency  resolution  of  the 
transmission  spectra  from  which  the  gaps  were  estimated,  (inset)  Magnitudes  of  the  same  gaps 
calculated  from  the  model  described  in  the  text. 
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Figure  3:  Optical  relaxation  gradients  at  avoided  crossings,  a  Measurements  of  (i)  cavity  relaxation 
rate  and  (ii)  detuning  for  each  cavity  mode  versus  position  for  the  membrane  tilted  0.66  mrad  about 
the  axis  indicated  in  the  figure,  b  Same  but  for  the  membrane  tilted  0.61  mrad  about  a  different 
axis  (indicated  in  the  figure).  The  dotted  arrows  highlight  the  crossing  whose  gradient  k/  reverses 
sign,  c  Relaxation  rate  k  of  the  TEMoo  mode  as  a  function  of  membrane  position.  Symbols  are 
data,  the  black  solid  line  is  a  smoothed  version  of  the  data  to  show  the  A/2  =  532  nm  periodicity. 
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Displacement  (nm) 


Figure  4:  Quartic  and  double-well  optomechanical  coupling,  a  Transmission  spectrum  for  the 
TEM20  mode  with  the  membrane  tilted  to:  (i)  1.0,  (ii)  1.3,  and  (iii)  1.6  mrad,  showing  transition 
from  (i)  quadratic  to  (ii)  quartic  to  (iii)  double- well  detunings,  b  Transmission  of  the  cavity  showing 
the  avoided  crossings  between  the  TEM{2o,ii,02}  and  TEM{4o73i722, 13,04}  modes. 
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Cooling  and  squeezing  via  quadratic  optomechanical  coupling 
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Departments  of  Physics  and  Applied  Physics,  Yale  University,  New  Haven,  Connecticut  06520,  USA 

(Dated:  April  21,  2010) 

We  explore  the  physics  of  optomechanical  systems  in  which  an  optical  cavity  mode  is  coupled  parametrically 
to  the  square  of  the  position  of  a  mechanical  oscillator.  We  derive  an  effective  master  equation  describing  two- 
phonon  cooling  of  the  mechanical  oscillator.  We  show  that  for  high  temperatures  and  weak  coupling,  the  steady- 
state  phonon  number  distribution  is  non-thermal  (Gaussian)  and  that  even  for  strong  cooling  the  mean  phonon 
number  remains  finite.  Moreover,  we  demonstrate  how  to  achieve  mechanical  squeezing  by  driving  the  cavity 
with  two  beams.  Finally,  we  calculate  the  optical  output  and  squeezing  spectra.  Implications  for  optomechanics 
experiments  with  the  membrane-in-the-middle  geometry  or  ultracold  atoms  in  optical  resonators  are  discussed. 

PACS  numbers:  37.30.+i,  42.50.Lc,  42.65. -k,  85.85.-hj 


Introduction.  In  optomechanical  systems  optical  and  me¬ 
chanical  degrees  of  freedom  are  coupled  via  radiation  pres¬ 
sure,  optical  gradient,  or  photo  thermal  forces.  While  work  in 
this  area  was  originally  motivated  by  the  goal  of  building  sen¬ 
sitive  detectors  for  gravitational  waves  [1,  2],  the  field  has  be¬ 
come  an  active  area  of  research  in  its  own  right.  Its  main  goal 
is  to  investigate  quantum  coherence  in  macroscopic  solid-state 
devices  both  for  quantum  information  purposes  and  gaining 
new  insights  into  the  quantum-to-classical  transition  [3,  4]. 

In  most  optomechanical  experiments  an  optical  cavity 
mode  is  parametrically  coupled  to  the  position  of  a  mechanical 
oscillator.  Consequently,  many  properties  of  this  setup  have 
been  discussed,  including  red-sideband  laser  cooling  in  the 
resolved- sideband  limit  [5,  6],  normal-mode  splitting  [7,  8], 
optical  squeezing  [9,  10],  backaction-evading  measurements 
[11-13],  mechanical  squeezing  using  either  feedback  [12, 14], 
squeezed  light  [15]  or  modulation  of  input  power  [16],  and  en¬ 
tanglement  between  light  and  a  mechanical  oscillator  [17]. 

However,  some  optomechanical  systems  feature  a  quadratic 
optomechanical  interaction,  i.e.  an  optical  cavity  mode  is  cou¬ 
pled  parametrically  to  the  square  of  the  position  of  a  mechan¬ 
ical  oscillator.  One  example  is  the  membrane-in-the-middle 
geometry,  in  which  the  membrane  is  placed  at  a  node  or  antin¬ 
ode  of  the  cavity  field  [18-20].  A  second  system  which  can 
realize  quadratic  optomechanical  coupling  is  a  cloud  of  ultra¬ 
cold  atoms  loaded  into  an  optical  cavity,  where  the  cloud’s 
center-of-mass  coordinate  serves  as  the  mechanical  degree  of 
freedom  [21].  To  date,  the  theoretical  literature  has  focused 
on  using  the  quadratic  coupling  to  detect  phonon  Fock  states 
[18,  19,  22-24],  but  otherwise  the  possible  uses  of  this  form 
of  optomechanical  coupling  are  largely  unstudied. 

In  this  paper  we  explore  three  features  of  quadratic  optome¬ 
chanical  coupling:  two-phonon  cooling  of  the  mechanical  os¬ 
cillator,  squeezing  of  the  mechanical  oscillator,  and  squeez¬ 
ing  of  the  optical  output  field.  Using  Fermi’s  Golden  rule  we 
first  write  down  an  effective  master  equation  for  the  mechan¬ 
ical  oscillator.  In  the  classical  limit  of  large  phonon  number, 
two-phonon  cooling  processes  change  the  steady-state  num¬ 
ber  distribution  from  exponential  to  Gaussian,  a  consequence 
of  the  nonlinear  damping.  In  the  quantum  limit  we  find  that 
ground-state  cooling  is  not  possible  since  two-phonon  cool¬ 


ing  processes  preserve  the  phonon-number  parity.  We  then 
demonstrate  that  the  model  maps  onto  a  degenerate  paramet¬ 
ric  oscillator  if  the  cavity  is  driven  by  two  laser  beams  whose 
frequencies  are  detuned  to  either  side  of  the  cavity  resonance 
by  an  amount  equal  to  the  mechanical  frequency.  This  opens 
up  the  possibility  of  mechanical  squeezing.  Finally,  we  cal¬ 
culate  the  optical  output  spectrum  and  find  that  this  system  is 
capable  of  producing  considerable  optical  squeezing. 

Hamiltonian.  We  start  from  the  Hamiltonian  (with  h  =  1) 

H  =  (ljr  +  gx 2)  (ala  -  (a^a))  +  ujM^b  +  H1  +  HK  (1) 

where  c or  is  the  cavity  resonance  frequency,  g  the  quadratic 
optomechanical  coupling,  and  x  =  xzpf(^  +  the  posi¬ 
tion  of  the  mechanical  oscillator  with  zero-point  fluctuations 
xzpf  =  (2mccM)-1^2,  frequency  ujm  and  mass  m.  a  and  b 
are  annihilation  operators  obeying  bosonic  commutation  re¬ 
lations.  H1  and  HK  describe  the  coupling  to  the  mechan¬ 
ical  and  optical  baths  and  the  optical  drive.  We  have  sub¬ 
tracted  the  steady-state  mean  photon  number  (aJa)  which 
renormalizes  the  frequency  of  the  mechanical  oscillator.  The 
Hamiltonian  (1)  is  relevant  to  systems  with  membrane-in- 
the-middle  geometry  [18],  to  ultracold  atoms  in  optical  res¬ 
onators  [21],  and  double-microdisk  whispering-gallery  mode 
resonators  [30]  when  the  first  derivative  of  the  cavity  disper¬ 
sion  relation  cccav(x)  vanishes,  i.e.  cc'av(xo)  =  0,  so  that 
g  =  cc''av(xo)/2  is  the  leading  order  of  the  optomechanical 
coupling. 

Expressing  a  =  e~luJLt{a  +  d)  with  the  laser  frequency  ul, 
choosing  a  real,  and  neglecting  (Fd  with  respect  to  a(dft  +  d), 
we  obtain  the  following  quantum  master  equation 


g  =  —i 


Hs,  Q 


+  KT>[d\g  +  7(1  +  nth)V[b\g  + ' inthV\b]]g 

(2) 


with  the  system  Hamiltonian 


Hs  =  -A  d)d  +  uMPb  +  g(b  +  S1)2^  +  d)  (3) 


where  A  =  u>l  —  ujr  is  the  detuning,  g  =  gax |PF  the  cou- 
pling,  k  and  7  the  cavity  and  mechanical  damping  rates,  and 
nt h  the  thermal  phonon  number.  We  assume  that  the  optical 
bath  is  at  zero  temperature.  V[o]g  =  6  go t  —  (o^og  +  go^o)/2 
denotes  the  standard  dissipator  in  Findblad  form. 


2 


FIG.  1:  Steady-state  mean  phonon  (Pb)  (left)  and  photon  number 
(dfd)  (right)  as  a  function  of  the  coupling  strength  g/n  and  thermal 
coupling  7/ft  for  thermal  phonon  number  nt h  =  1,  obtained  from 
the  numerical  solution  of  the  full  quantum  master  equation  (2). 


Numerical  simulations  of  the  full  quantum  master  equa¬ 
tion.  If  the  cavity  is  driven  on  the  red  two-phonon  resonance, 
i.e.  A  =  —2ujm,  in  the  good-cavity  limit,  i.e.  ft  <C  c om,  we 
expect  two-phonon  cooling  processes  to  be  important.  Con¬ 
centrating  on  the  resonant  terms,  i.e.  Hs  =  g(PPd  +  h.c.), 
we  solve  the  full  quantum  master  equation  (2)  numerically  for 
small  thermal  phonon  numbers  nth-  We  find  three  different 
regimes:  (i)  for  7nth  ft  the  mechanical  oscillator  is  cooled 
due  to  the  coupling  to  the  zero-temperature  bath  of  the  optical 
field  [25],  (ii)  for  7nth  +  ft  the  optical  field  is  heated  by  the 
coupling  to  the  mechanical  oscillator,  and  (iii)  for  7nth  ~  ft 
both  effects  are  important  and  the  density  matrix  has  non-zero 
off-diagonal  elements  so  that  e.g.  the  correlator  \(cftbb)\  is 
non-zero.  As  an  example  we  plot  in  Fig.  1  the  steady- state 
mean  phonon  number  (m  b)  and  the  mean  number  of  photons 
due  to  the  coupling  to  the  membrane  (wd)  as  a  function  of  the 
coupling  strength  g/n  and  thermal  coupling  7/ ft  for  thermal 

phonon  number  7ith  =  1. 

Effective  master  equation  describing  two-phonon  cooling. 
For  771th  ft  and  weak  coupling,  i.e.  g  <C  ft,  we  can  em¬ 
ploy  a  quantum  noise  approach  [25],  i.e.  we  calculate  two- 
phonon  cooling  and  amplification  rates  using  Fermi’s  Golden 
rule.  Concentrating  on  the  diagonal  terms  of  the  density  ma¬ 
trix  gnn  =  Pn,  we  write  down  a  set  of  rate  equations 

P n  =  —7  (ftthfft  +  1)  +  (ftth  +  l)ft)  Pn 
+  'ynthnPn-i  +  7(nth  +  l)(n  +  l)Pn+i 
—  (F+(n  —  1)  +  rt(n  +  2)(n  + 1))  pn 

+  I\(ft  +  2  )(n  +  l)Pn+2  +  r^n(n  —  l)Pn_2*  (4) 

The  terms  in  the  first  two  lines  are  due  the  coupling  of  the  me¬ 
chanical  oscillator  to  its  thermal  bath  with  rate  7  and  thermal 
phonon  number  nth-  The  terms  in  the  last  two  lines  are  due 
to  two-phonon  processes  whose  rates  are  governed  by  = 
g2x^¥FSnn(±2uM)  where  Snn(w)  =  k,\cl\2\xr{u)\2  is  the 
photon  number  spectral  density  and  xr( w)  =  [ft/ 2  —  i(cc  + 
A)]-1  the  cavity  response  function.  We  note  that  nonlinear 
damping  of  a  mechanical  oscillator  was  studied  in  Ref.  [26] . 

The  infinite  set  of  rate  equations  (4)  can  be  replaced 
by  a  single  differential  equation  for  the  generating  function 
F(z,  t )  =  Pn{t)zn-  For  vanishing  two-phonon  ampli¬ 

fication  Tt  =  0  the  steady- state  equation  can  be  solved  ex¬ 
actly  in  terms  of  the  confluent  hypergeometric  function  [27] . 


FIG.  2:  Steady- state  phonon  number  distribution  Pn  for  thermal 
phonon  number  nt h  =  107  and  for  =  0  (dashed)  and  r^/7  = 
4  •  10-7  (solid),  respectively. 


We  point  out  that  it  is  an  advantage  of  quadratic  cooling  that 
it  enables  cooling  when  the  membrane  or  atoms  are  placed  at 
a  node  of  the  cavity  field.  In  this  situation  the  system  is  most 
insensitive  to  absorption  from  the  membrane  or  atoms  as  well 
as  to  bistability  from  radiation  pressure. 

Two-phonon  cooling  in  the  classical  limit.  For  large  ther¬ 
mal  phonon  number  nth  +  1  and  weak  two-phonon  cooling 
771th  +  T;  we  can  replace  the  quantum  operators  b  and  d 
in  their  Heisenberg  equations  of  motion  by  complex  ampli¬ 
tudes  f3  =  (b)  and  a  =  (d)  and  obtain  two  coupled  clas¬ 
sical  Langevin  equations  $  =  —7/3/2  —  2 igf3*a  +  £  and 
a  =  —na/2  —  ig/3 2,  where  the  thermal  noise  is  characterized 
by  (O)  =  7nth S(t  —  t').  Adiabatically  eliminating  the 

optical  field  we  obtain  an  equation  of  motion  with  a  nonlinear 
damping  term  /3  =  —7/3/2  —  Ag2\f3\2 (3 / n  +  £.  Solving  the 
corresponding  Fokker-Planck  equation  we  obtain  the  phonon 
number  distribution 


Pn  oc  exp 


(5) 


The  distribution  changes  from  an  exponential  for  7/r^  nth 
to  a  Gaussian  for  7/^  <C  ftth-  In  the  latter  limit  the  mean 
phonon  number  is  given  by  (n)  =  y/7ftthft/7r^2-  The  ex¬ 
pression  in  Eq.  (5)  agrees  with  the  high-temperature  limit  of 
the  exact  solution  to  the  rate  equations  [27] .  We  conclude  that 
the  change  in  the  steady- state  phonon  number  distribution  is  a 
purely  classical  effect  due  to  nonlinear  damping. 

The  quadratic  coupling  in  current  membrane-in-the-middle 
experiments  is  very  small  g/n  =  10-5.  Nonetheless,  it  leads 
to  sizable  effects  if  the  thermal  phonon  number  nt h  is  large. 
In  Fig.  2  we  plot  the  steady- state  phonon  number  distribution 
for  a  thermal  phonon  number  nth  =  107  both  in  the  absence 
T^  =  0  and  presence  r^/7  =  4  •  10“ 7  of  two-phonon  cooling. 
We  find  the  Planck  distribution  with  mean  nt h  =  107  becomes 
a  nearly-Gaussian  distribution  with  mean  nt h  =  2.8  •  106. 

Two-phonon  cooling  in  the  quantum  limit.  Experiments 
with  ultracold  atoms  are  in  the  opposite  limit  with  small  ther¬ 
mal  heating  rate  771th  and  strong  optomechanical  coupling  g. 
Solving  the  rate  equations  (4)  in  the  limit  of  strong  optical 
damping,  i.e.  T^  777th ,  we  obtain  the  phonon  number  dis¬ 
tribution  Pq/Pi  =  3  +  1/  77th  with  all  other  Pn  =  0  and  the 
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FIG.  3:  (Color  online)  Mean  phonon  number  (Pb)  as  a  function  of 
coupling  g/u  obtained  from  the  quantum  master  equation  (blue  dots) 
and  the  classical  rate  equations  without  (red  solid)  and  with  strong¬ 
coupling  correction  (green  dashed).  The  black  dashed  line  indicates 
the  minimal  phonon  number  n°M  given  by  Eq.  (6)  for  nt h  =  5. 


minimal  mean  phonon  number 


fjO  _ 
nM  ~ 


4  -f — — 

nth 


(6) 


The  fact  that  strong  cooling  leaves  both  ground  and  first  ex¬ 
cited  state  occupied  is  a  consequence  of  the  fact  that  two- 
phonon  cooling  processes  preserve  the  phonon-number  parity. 

In  Fig.  3  we  plot  the  steady-state  phonon  number  (mb)  as  a 
function  of  coupling  g  obtained  from  the  full  quantum  master 
equation  (2)  and  the  rate  equations  (4).  The  two  results  show 
excellent  agreement  at  weak  coupling  g  <C  ft.  As  4g2/ft  be¬ 
comes  comparable  to  ft,  the  Fermi’s  Golden  rule  expression 
for  two-phonon  cooling  breaks  down  and  the  predictions  of 
the  full  master  equation  (2)  and  the  weak-coupling  rate  equa¬ 
tions  (4)  will  in  general  be  different.  In  the  limit  where  only 
the  Fock  states  with  phonon  number  n  <  3  are  important, 
we  adiabatically  eliminate  the  off-diagonal  terms  in  the  quan¬ 
tum  master  equation  (2)  and  find  that  the  strong-coupling  two- 
phonon  cooling  rates  are  given  by 


■  n,n— 2 


4  g2n(n  —  l)ft 
ft2  +  4  g2n(n  —  1) 


(7) 


In  the  weak-coupling  limit  g  <C  ft  this  expression  simplifies 
to  our  previous  result  4 g2n(n  —  l)/ft.  In  the  limit  of  strong 
coupling  g  ft,  the  two-phonon  cooling  rate  n_2  remains 
finite  and  cannot  exceed  the  cavity  damping  rate  ft.  This  leads 
to  a  minimum  phonon  number  which  is  larger  than  the  one 
predicted  from  the  weak-coupling  theory.  In  Fig.  3  we  plot  the 
steady- state  mean  phonon  number  obtained  from  the  modified 
rate  equations  and  find  excellent  agreement  with  the  exact  so¬ 
lution  to  the  full  quantum  master  equation  (2). 

We  note  that  for  quadratic  coupling  in  the  good-cavity  limit, 
the  phonon  number  distribution  Pn  can  be  measured  by  mon¬ 
itoring  the  phase  shift  of  the  reflected  light  [18,  19,  22]. 

Squeezing  the  mechanical  oscillator.  Driving  the  cavity  at 
both  ujr  ±  ujm  with  equal  strength,  the  classical  part  of  the 
cavity  field  oscillates  in  time  a  =  A  cos  In  the  case 


FIG.  4:  Optimal  squeezing  spectrum  of  the  cavity  output  field 
S'opt  (cj)  as  a  function  of  detuning  A.  The  parameters  are  uJm/n  =  5, 
g/u  —  0.5,  7/ft  —  0.1  and  nth  —  0. 


of  linear  optomechanical  coupling  this  enables  a  backaction- 
evading  measurement  of  one  quadrature  of  the  mechanical  os¬ 
cillator  [11-13],  but  does  not  itself  produce  squeezing  apart 
from  the  one  which  is  conditioned  on  the  measurement  out¬ 
come.  When  we  instead  consider  the  same  drive  applied  to 
a  system  with  quadratic  coupling,  moving  to  an  interaction 
picture  with  respect  to  Hq  =  urcl^cl  +  utM^b  and  keep¬ 
ing  only  non-rotating  terms  of  the  quadratic  coupling,  we  ob¬ 
tain  the  Hamiltonian  of  the  degenerate  parametric  oscillator, 
i.e.  Hd po  =  f  (b2  +  (Sf)2)  with  x  =  gx\PFA/2. 

Solving  the  linear  quantum  Langevin  equations  we  see  that 
the  steady- state  fluctuations  in  the  quadrature  X  =  (beluJMt  + 
H.c.)/  y/2  of  the  mechanical  oscillator  are  squeezed  below  the 
thermal  level,  (X2)  =  ’  depending  on  the  ratio  x/l- 

At  threshold  x  =  7/2,  i.e.  before  the  parametric  oscillator 
becomes  unstable,  it  is  maximal  and  equal  to  -3dB  [28].  Al¬ 
though  the  parameters  of  current  membrane-in-the-middle  se¬ 
tups  show  small  coupling  g  <C  ft  and  large  thermal  phonon 
number  nt h  1,  we  emphasize  that  the  ratio  x/l  which  is 

of  importance  here  can  still  be  comparable  to  unity  and  lead 
to  significant  noise  squashing. 

The  mechanical  squeezing  can  be  detected  by  coupling  the 
position  of  the  mechanical  oscillator  parametrically  to  a  sec¬ 
ond  optical  mode.  This  scenario  has  been  studied  in  Ref.  [14]. 

Output  spectrum.  Let  us  now  return  to  the  full  model  (3). 
Up  to  second  order  in  the  coupling  g ,  the  cavity  output  spec¬ 
trum  which  is  defined  as  (cc)  =  f  dt  eluJt  (dQ ut  (t)d/out  (0)) 
with  dout  =  din  +  \Z~dd  and  (d[n(t)dln(t))  =  S(t  - 1')  is  given 
by  sdd(w)  =  4k92\Xr[-u]\2Sx2x2(u)  where 


& X2X 2  7) 


7(nth  + 1)2 

7 2  +  (w  —  2  wm)2 

7»th _ 27nth(nth  +  1) 

72  +  (w  +  2  lum)2  72  +  w2 


(8) 


We  see  the  output  spectrum  (w)  has  sidebands  at  oj  = 
±2 Lum  and  uj  =  0  as  expected  for  quadratic  coupling. 

Optical  squeezing  spectrum.  One  application  of  optome¬ 
chanical  devices  which  has  been  widely  advocated  [9,  10] 
is  to  use  the  nonlinear  coupling  between  light  and  mirror  to 
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squeeze  the  incoming  coherent  light  beam,  i.e.  reduce  one  of 
its  quadratures  below  the  shot-noise  level  at  certain  frequen¬ 
cies.  The  quantity  which  characterizes  this  noise  reduction  is 
the  optical  squeezing  spectrum  SffUt(cu)  given  by  [29] 


tM  =  1  +  J 


dtelut{:  X%u\t),X%ut( 0)  :) 


—  OO 
OO 


=  1  +  k  f  dt  e*wtT 

J  —  OO 


(:Xe(t),Xe( 0)  :) 


(9) 


where  Xgut  =  (Souteie+H.c.)/ 2,  (A,  B)  =  (AB)-(A)(B), 

the  colons  indicate  normal  ordering  and  T  is  the  time-ordering 
operator  [29] .  The  former  expression  is  useful  for  calculations 
in  input-output  theory  and  the  latter  for  master  equation  sim¬ 
ulations  evoking  the  quantum  regression  theorem. 

In  Fig.  4  we  plot  the  optimal  squeezing  spectrum  of  the  cav¬ 
ity  output  field  (cu)  =  min#  S^ut(cu)  as  a  function  of  de¬ 
tuning  A.  In  the  good-cavity  limit  at  zero  temperature  we  find 
considerable  squeezing  at  uj  =  0  and  uj  =  2ujm  for  A  =  0  and 
A  =  —  2ccm>  respectively.  We  choose  parameters  ccm/ft  =  5, 
g/n  =  0.5,  7 j k  =  0.1  and  nth  =  0,  relevant  to  experiments 
with  ultracold  atoms  in  optical  resonators. 

To  gain  insight  beyond  numerics  we  obtain  the  squeezing 
spectrum  perturbatively  up  to  second  order  in  the  coupling  g 


s™\(u)  =  i  +  sz\u)  +  sz\-u) 

-  8g2 k,\xr(u)xr(~u)Sx2x2(u) 


—  ^th  + 


X*r(u)xr(~u)  (kXaM  -  1) 

1 


7  +  i(cu  -  2u>m)  7  +  i(oj  +  2u>m) 


•  (10) 


As  in  the  case  of  linear  optomechanical  coupling  [9,  10],  we 
find  two  regimes  of  strong  squeezing:  for  small  detuning  and 
small  frequencies  as  well  as  for  detuning  and  frequencies 
close  to  twice  the  mechanical  frequency.  However,  thermal 
fluctuations  are  more  destructive  for  squeezing  effects  with 
quadratic  optomechanical  coupling  due  to  the  quadratic  scal¬ 
ing  with  the  thermal  phonon  number  nth  and  due  to  the  fact 
that  the  spectrum  SX2X2  (uj)  has  weight  at  zero  frequency. 

Conclusions.  We  have  explored  the  physics  of  nonlinear 
optomechanical  systems  where  an  optical  cavity  mode  cou¬ 
ples  quadratically  rather  than  linearly  to  the  position  of  a  me¬ 
chanical  oscillator.  For  optomechanical  experiments  with  the 
membrane-in-the-middle  geometry,  we  predict  a  qualitative 
change  in  the  phonon  number  distribution  and  mechanical 
noise  squashing.  For  ultracold  atoms  in  optical  resonators,  we 
found  the  quantum  limit  for  two-phonon  cooling  and  sizable 
mechanical  and  optical  squeezing. 
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